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An interesting line of research is the investigation of the laws of random variables known as 
Dirichlet means. However, there is not much information on interrelationships between different 
Dirichlet means. Here, we introduce two distributional operations, one of which consists of 
multiplying a mean functional by an independent beta random variable, the other being an 
operation involving an exponential change of measure. These operations identify relationships 
between different means and their densities. This allows one to use the often considerable analytic 
work on obtaining results for one Dirichlet mean to obtain results for an entire family of otherwise 
seemingly unrelated Dirichlet means. Additionally, it allows one to obtain explicit densities for 
the related class of random variables that have generalized gamma convolution distributions 
and the finite-dimensional distribution of their associated Levy processes. The importance of 
this latter statement is that Levy processes now commonly appear in a variety of applications 
in probability and statistics, but there are relatively few cases where the relevant densities have 
been described explicitly. We demonstrate how the technique allows one to obtain the finite- 
dimensional distribution of several interesting subordinators which have recently appeared in 
the literature. 

Keywords: beta-gamma algebra; Dirichlet means and processes; exponential tilting; 
generalized gamma convolutions; Levy processes 

1. Introduction 

In this work, we present two distributional operations which identify relationships be- 
tween seemingly different classes of random variables which are representable as linear 
functionals of a Dirichlet process, otherwise known as Dirichlet means. Specifically, the 
first operation consists of multiplication of a Dirichlet mean by an independent beta ran- 
dom variable and the second operation involves an exponential change of measure to the 
density of a related infinitely divisible random variable having a generalized gamma con- 
volution distribution (GGC). This latter operation is often referred to in the statistical 
literature as exponential tilting or in mathematical finance as an Esscher transform. We 
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believe our results add a significant component to the foundational work of Cifarelli and 
Regazzini [5, 6]. In particular, our results allow one to use the often considerable analytic 
work on obtaining results for one Dirichlet mean to obtain results for an entire family 
of otherwise seemingly unrelated mean functionals. It also allows one to obtain explicit 
densities for the related class of infinitely divisible random variables which are generalized 
gamma convolutions and an explicit description of the finite-dimensional distribution of 
their associated Levy processes (see Bertoin [1] for the formalities of general Levy pro- 
cesses). The importance of this latter statement is that Levy processes now commonly 
appear in a variety of applications in probability and statistics, but there are relatively 
few cases where the relevant densities have been described explicitly. A detailed summary 
and outline of our results may be found in Section 1.2. Some background information on, 
and notation for, Dirichlet processes and Dirichlet means, their connection with GGC 
random variables, recent references and some motivation for our work are given in the 
next section. 

1.1. Background and motivation 

Let X be a non- negative random variable with cumulative distribution function Fx- 
Note, furthermore, that for a measurable set C, we use the notation Fx{C) to mean 
the probability that X is in C. One may define a Dirichlet process random probabil- 
ity measure (see Freedman [17] and Ferguson [15, 16]), say Pg, on [0,oo) with total 
mass parameter 9 and prior parameter Fx, via its finite-dimensional distribution as 
follows: for any disjoint partition on [0, oo), say (C\, . . . , Ck), the distribution of the ran- 
dom vector (Pg(Ci), . . . ,Pg(Ck)) is a fc-variate Dirichlet distribution with parameters 
(0F X (G 1 ), ...,6F x (C k )). Hence, for each C, 



has a beta distribution with parameters (9 Fx (C), 0(1 — Fx(C))). Equivalently, setting 
6F x (C i ) = 6 i for i = l,..., k, 



where (GgJ are independent random variables with gamma(#i, 1) distributions and Gg = 
Gg 1 + ■■■ + Gg k has a gamma(0, 1) distribution. This means that one can define the 
Dirichlet process via the normalization of an independent increment gamma process on 
[0, oo), say ^g(-), as 





7fl([0,oo)) 



where 7e(C,) = G$ i} whose almost surely finite total random mass is 7e([0,oo)) = Gg. A 
very important aspect of this construction is the fact that Gg is independent of Pg and 
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hence of any functional of Pq. This is a natural generalization of Lukacs' [35] characteriza- 
tion of beta and gamma random variables, which is fundamental to what is now referred 
to as the beta-gamma algebra (for more on this, see Chaumont and Yor ([4], Section 4.2); 
see also Emery and Yor [11] for some interesting relationships between gamma processes, 
Dirichlet processes and Brownian bridges). Hereafter, for a random probability measure 
P on [0,oo), we write 



to indicate that P is a Dirichlet process with parameters (9, Fx)- 

These simple representations and other nice features of the Dirichlet process have, 
since the important work of Ferguson [15, 16], contributed greatly to the relevance and 
practical utility of the field of Bayesian non- and semi-parametric statistics. Naturally, 
owing to the ubiquity of the gamma and beta random variables, the Dirichlet process 
also arises in other areas. One of the more interesting and, we believe, quite important 
topics related to the Dirichlet process is the study of the laws of random variables called 
Dirichlet mean Junctionals, or simply Dirichlet means, which we denote as 



as initiated in the works of Cifarclli and Regazzini [5, 6]. In [6], the authors obtained an 
important identity for the Cauchy-Stieltjes transform of order 9. This identity is often 
referred to as the Markov-Krein identity, as can be seen in, for example, Diaconis and 
Kcmperman [9], Kerov [28] and Vershik, Yor and Tsilevich [40], where these authors 
highlight its importance to, for instance, the study of the Markov moment problem, con- 
tinued fraction theory and exponential representation of analytic functions. This identity 
is later called the Cifarelli-Regazzini identity in [21]. Cifarclli and Regazzini [6], owing to 
their primary interest, used this identity to then obtain explicit density and cdf formulae 
for Mg (Fx ) ■ The density formulae may be seen as Abel- type transforms and hence do not 
always have simple forms, although we stress that they are still useful for some analytic 
calculations. The general exception is the case = 1, which has a nice form. Some exam- 
ples of works that have proceeded along these lines arc Cifarclli and Mclilli [7], Regazzini, 
Guglielmi and di Nunno [38], Regazzini, Lijoi and Priinster [39], Hjort and Ongaro [20], 
Lijoi and Regazzini [32], and Epifani, Guglielmi and Melilli [12, 13]. Moreover, the recent 
works of Bertoin et al. [2] and James, Lijoi and Priinster [25] (see also [23], which is 
a preliminary version of this work) show that the study of mean functionals is relevant to 
the analysis of phenomena related to Besscl and Brownian processes. In fact, the work of 
James, Lijoi and Priinster [25] identifies many new explicit examples of Dirichlet means 
which have interesting interpretations. 

Related to these last points, Lijoi and Regazzini [32] have highlighted a close connec- 
tion to the theory of generalized gamma convolutions (see [3]). Specifically, it is known 
that a rich subclass of random variables having generalized gamma convolutions (GGC) 
distributions may be represented as 





(1.1) 
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We call these random variables GGC(6*,Fx). In addition, we see from (1.1) that Tg 
is a random variable derived from a weighted gamma process and, hence, the calculus 
discussed in Lo [33] and Lo and Wcng [34] applies. In general, GGC random variables are 
an important class of infinitely divisible random variables whose properties have been 
extensively studied by [3] and others. We note further that although we have written a 
GGC(#, Fx) random variable as GgMg(Fx), this representation is not unique and, in 
fact, it is quite rare to see Tg represented in this way. We will show that one can, in fact, 
exploit this non- uniqueness to obtain explicit densities for Tg , even when it is not so easy 
to do so for Mg(Fx)- While the representation GgMg(Fx) is not unique, it helps one 
to understand the relationship between the Laplace transform of Tg and the Cauchy- 
Stieltjes transform of order 8 of Mg(Fx), which, indeed, characterize respectively the laws 
of Tg and Mg(Fx)- Specifically, using the independence property of Gg and Mg(Fx) leads 
to, for A > 0, 

E[e~ ATe ] =E[(l + AM„(F x ))^]=c^ W , (1-2) 

where 

poo 

^f x (X)= \og(l + \x)F x (dx)=E[\og(l + \X)} (1.3) 
Jo 

is the Levy exponent of Tg. We note that Tg and Mg(Fx) exist if and only if ipp x (A) < oo 
for A > (see, e.g., [8] and [3]). The expressions in (1.2) equate with the aforementioned 
identity obtained by Cifarelli and Regazzini [6]. 

Despite these interesting results, there is very little work on the relationship between 
different mean functionals. Suppose, for instance, that for each fixed value of 6 > 0, 
Mg(Fx) denotes a Dirichlet mean and (Mg(Fz c );c > 0) denotes a collection of Dirichlet 
mean random variables indexed by a family of distributions (Fz e ; c> 0). One may then 
ask the following question: for what choices of X and Z c are these mean functionals 
related, and in what sense? In particular, one may wish to know how their densities are 
related. The rationale here is that if such a relationship is established, then the effort that 
one puts forth to obtain results such as the explicit density of Mg(Fx) can be applied to 
an entire family of Dirichlet means (Mg(Fz c ); c> 0). Furthermore, since Dirichlet means 
are associated with GGC random variables, this would establish relationships between 
a GGC(6>, Fx) random variable and a family of GGG(6,Fz c ) random variables. Simple 
examples are, of course, the choices Z c = X + c and Z c = cX, which, due to the linearity 
properties of mean functionals, result easily in the identities in law 

Mg{F x+ c) = c + M e (Fx) and Mg (F cX ) = cMg {F x ) . 

Naturally, we are going to discuss more complex relationships, but with the same goal. 
That is, we will identify non-trivial relationships so that the often considerable efforts that 
one makes in the study of one mean functional Mg (Fx ) can then be used to obtain more 
easily results for other mean functionals, their corresponding GGC random variables and 
Levy processes. In this paper, we will describe two such operations which we elaborate 
on in the next subsection. 
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1.2. Outline and summary of results 

Section 1.3 reviews some of the existing formulae for the densities and cdfs of Dirichlet 
means. In Section 2, we will describe the operation of multiplying a mean functional 
Mda(Fx) by an independent beta random variable with parameters (9a, 9(1 — a)), say, 
Pe<7,e(i-cr)i where < a < 1. We call this operation beta scaling. Theorem 2.1 shows 
that the resulting random variable f3g a ^g(i- a )Mg a (Fx) is again a mean functional, but 
now of order 9. In addition, the GGC(9<j, Fx) random variable Gg a Mg a (Fx) is equiv- 
alcntly a GGC random variable of order 9. Now, keeping in mind that tractable den- 
sities of mean functionals of order 9 = 1 are the easiest to obtain, Theorem 2.1 shows 
that by setting 9 = 1, the densities of the uncountable collection of random variables 
(f3 (J: i^ cr M a (Fx)', < a < 1) are all mean functionals of order 9 = 1. Theorem 2.2 then 
shows that efforts used to calculate the explicit density of any one of these random vari- 
ables, via the formulae of [6] , lead to the explicit calculation of the densities of all of them. 
Additionally, Theorem 2.2 shows that the corresponding GGC random variables may all 
be expressed as GGC random variables of order 9 = 1, representable in distribution as 
G\fi a ^- ( yM ( ,(Fx)- A key point here is that Theorem 2.2 gives a tractable density for 
{i<j,\- a M a (Fx) without requiring knowledge of the density of M„(Fx), which is usually 
expressed in a complicated manner. These results will also yield some non-obvious inte- 
gral identities. Furthermore, noting that a GGC (9, Fx) random variable, Tg, is infinitely 
divisible, we associate it with an independent increment process (£g(t);t > 0), known as 
a subordinator (a non-decreasing non- negative Levy process), where, for each fixed t, 

E[e- A ^ (t >] =E[c- AT ^] = c -Wr x M. 

That is, marginally, Ce(l) = Tg and (g(t) = Cet(l) = Tg t . In addition, for s < t, Ce(t) — 
(g(s) = (g(t — s) is independent of (g(s). We say that the process ((g(t);t > 0) is a 
GGC(9, Fx) subordinator. Proposition 2.1 shows how Theorems 2.1 and 2.2 can be 
used to address the usually difficult problem of explicitly describing the densities of 
the finite-dimensional distribution of a subordinator (see [29]). This has implications in, 
for instance, the explicit description of densities of Bayesian nonparametric prior and 
posterior models, but is clearly of wider interest in terms of the distribution theory of 
infinitely divisible random variables and associated processes. 

In Section 3, we describe how the operation of exponentially tilting the density of 
a GGC (9, Fx) random variable leads to a relationship between the densities of the mean 
functional Mg(Fx) and yet another family of mean functionals. This is summarized in 
Theorem 3.1. Section 3.1 then discusses a combination of the two operations. Propo- 
sition 3.1 describes the density of beta-scaled and tilted mean functionals of order 1. 
Using this, Proposition 3.2 describes a method to calculate a key quantity in the explicit 
description of the densities and cdfs of mean functionals. In Section 4, we show how 
the results in Sections 2 and 3 are used to derive the finite-dimensional distribution and 
related quantities for classes of subordinators suggested by the recent work of James, 
Lijoi and Prunster [25] and Bertoin et al. [2]. 
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1.3. Preliminaries 



Suppose that X is a positive random variable with distribution Fx and define the function 

/•OO 

/ log(|t-x|)I(Mx)F x (d a ;)=E[lo g (|t-X|)I(<^X)]. 



Furthermore, define 

A e (t|F x ) = isin(7t^xW)c" e * Fx( * ) , 

where, using a Lebesgue-Stieltjes integral, Fx(t) = f Q Fx{dx). Cifarelli and Regazz- 
ini [6] (see also [7]) apply an inversion formula to obtain the distributional formula for 
Mq(Fx) as follows. For all 9 > 0, the cdf can be expressed as 

(x-i) e_1 A e (i|F x )dt, (1.4) 

provided that 9 Fx possesses no jumps of size greater than or equal to one. If we let 
£,8F X (') denote the density of Mg(Fx), then it takes its simplest form for 9 = 1, which is 

t Fx (x) = ArixlFx) = -sm(nFx(x))c-^. (1.5) 

Density formulae for 9 > 1 are described as 

^ BFx {x) = {6-l) ( (x-t) 9 ~ 2 Ae(t\F x )dt. (1.6) 



An expression for the density, which holds for all 9 > 0, was recently obtained by 
James, Lijoi and Priinstcr [25] as follows: 

Z8F X {x) = - f (x-tf^deitlFx)^, (1.7) 
n Jo 

where 

d e (t\F x ) = ± S m(n6Fx{t))e- e *W. 
For additional formulae, see [6, 32, 38]. 

Remark 1.1. Throughout, for random variables R and X, when we write the product 
RX, we will assume, unless otherwise mentioned, that R and X are independent. This 
convention will also apply to the multiplication of the special random variables that are 
expressed as mean functionals. That is, the product Mg(Fx)Mg(Fz) is understood to be 
a product of independent Dirichlet means. 
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Remark 1.2. Throughout, we will be using the fact that if R is a gamma random 

variable, then the independent random variables R, X, Z satisfying RX = RZ imply that 

X = Z. This is true because gamma random variables arc simplifiablc. For the precise 
meaning of this term and associated conditions, see Chaumont and Yor [4], Sections 1.12 
and 1.13. This fact also applies to the case where R is a positive stable random variable. 

2. Beta scaling 

In this section, we investigate the simple operation of multiplying a Dirichlet mean 
functional Mg(Fx) by certain beta random variables. Note, first, that if M denotes an 
arbitrary positive random variable with density /m, then, by elementary arguments, it 
follows that the random variable W = fi a ^M, where is beta(a, b) independent of M, 
has density expressible as 

= Tr^rk f fM(w/u)u a - 2 (i uf- 1 du. 

T(a)T(b) Jo 

However, it is only in special cases that the density fw can be expressed in even simpler 
terms. That is to say, it is not obvious how to carry out the integration. In the next 
results, we show how remarkable simplifications can be achieved when M = Mg(Fx), in 
particular, for the range < 6 < 1, and when j3 a ^ is a symmetric beta random variable. 
First, we will need to introduce some additional notation. Let Y a denote a Bernoulli 
random variable with success probability < a < 1. Then, if X is a random variable 
with distribution Fx , independent of Y a , it follows that the random variable XY a has 
distribution 

F X yMx) = vF x (dx) + (1 - a)d (dx) (2.1) 

and cdf 

F XYa (x) = aFx(x) + (1 - a)I(x > 0). (2.2) 
Hence, there exists the mean functional 

M e (F XYa )= / yPe(dy), 
Jo 

where Pg(dy) denotes a Dirichlet process with parameters (0,Fxy„)- In addition, we 
have, for x > 0, 

•$>F XYa 0) = Epogfli - XY a \)l{XY a ± x)} = a$ Fx (x) + (1 - a) log(x). (2.3) 

When a = 1, Y a = 1 and hence FxyA?) = Fx{x). Let E a denote a set such that 
¥\Pg(E a )] = a. Note, now, that every beta random variable, /3 a .b, where a,b are arbi- 
trary positive constants, can be represented as the simple mean functional 



Pe(E a )±p ecA i-„) = Me{FY„), 
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by choosing 

a 

a = and 9 = a + b. 

a + b 

We note, however, that there are other choices of Fx that will also yield beta random 
variables as mean functionals. Throughout, we will use the convention that Pe,o '■= lj 
that is, the case where a = 1. We now present our first result. 

Theorem 2.1. For 9 > and < a < 1, Zei Pg a .e(\-a) denote a beta random variable 
with parameters (9a, 9(1 — cr)), independent of the mean functional Mg a (Fx)- Then: 

(i) Pe„,e{i-<,)M 6< r(F x ) = Mg(F X Y„)\ 

(ii) equivalents, Mg (F Ya )Mg a (F x ) = Mg (F X Y a ) \ 

(iii) Gg a Mg a (F x ) = G e Mg(F XYa ); 

(iv) that is, GGC(9a 1 Fx)=GGC(9 1 F XY „). 

Proof. Since Me(Fy CJ ) = Peafi{\-o)i statements (i) and (ii) are equivalent. We proceed 
by first establishing (iii) and (iv). Note that, using (1.3), 

/•OO 

E[log(l + XXY a )] = CT E[log(l + XX)} = cr / log(l + Xx)F x (dx). 

Jo 

Hence, 

^r e -XG e M e (F XYa )-i _ c -0<y f~ \o g (l+\x)F x (dx) = jg^-AG^M^ (F x ) j ^ 

which means that GgMg(FxY CJ ) = Gig a Mg a (Fx), establishing statements (iii) and (iv). 
Now, writing Gg a = Ggf3g a ^(i_,j), h follows that 

GgMg(FxY„) = Gg/3g a ^^i_ cr ^Mg a (Fx)- 

Hence, Pg a .e(i-a)Mg a (Fx) = Mg(FxY a ), by the fact that gamma random variables are 
simplifiablc. □ 

Remark 2.1. Wc note that parts (i) and (ii) of Theorem 2.1 also follow as consequences 
of Ethier and Griffiths [14], Lemma 1. We now state their interesting result for clarity. 

Lemma 2.1 (Ethier and Griffiths [14]). Let v\ and denote two probability mea- 
sures. Now, for #i,f?2 > 0, define the probability measure 

9 9 

V( 01 ,g 2 )(dx) = Q ^ V!(Ax) + ^ V2(dx). 

Then, for independent Dirichlet processes fj,± ~ Iig 1 ^ 1 and /12 ~ Hg 2tU2 , 
Mi,2(0 = Pe u e 2 (J-i(-) + (1 - /3ei,02 W'): 
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where fii^2 is a Dirichlet process with parameters (6\ + 82, Vfe^fo))- 
Hence, as a general consequence, 

Me 1+ e 2 (v(e u e 2 )) = Pe u e 2 M ei{vi) + (1 - /3e 1: e 2 )Me 2 {v2)- 

Now, from (2.1), we see that setting V\ = Fx, 1^2 = So, 61 = da and 62 = 0(1 — a) yields 
statements (i) and (ii). This is because Me{\- a )(So) = 0. 

When 6 — 1, wc obtain results for random variables (3 G ^\- a M a (Fx)- The symmetric 
beta random variables fia,\—a arise in a variety of important contexts and are often 
referred to as generalized arcsine laws with density expressible as 

*j^u°-\l-uy° for0<u<l. 
Now, using (2.1) and (2.2), let C{F X ) = {x : F x (x) > 0}. Then, for x > 0, 

sin(nF X Y (x)) = { Shi (™l 1 -^ x )V> ?**%i X \' (2.4) 
v aX " \sm(7t(l -a)), \ix£C(F x ). 

Also, note that sin(7t[l — F x (x)]) = sm(nF x (x)). The next result yields another surprising 
property of these random variables. 



Theorem 2.2. Consider the setting in Theorem 2.1. Then, when 8~ 1, it follows that 
for each fixed < a < 1, the random variable Mi(F X Y a ) = Pa.i-aMa{F x ) has density 

^ Fx ^(x) = ^ r sm(nF XYcr (x))e-^^ forx>0, (2.5) 

specified by (2.4)- Since GGC(o~, Fx) = GGC(l,FxY a ), this implies that the random vari- 
able G a M a (F x ) = G\M 1 (Fxy„) has density 



1 f°° 

g a , Fx (x) = - / e"*/ V" 2 sm(nF XYa (y))e" ff *-x M dy. 



(2.6) 



Proof. Since Mi(FxY cr ) = Pa.i-aMa{F x ), the density is of the form (1.5) for each fixed 
a e (0, 1]. Furthermore, we use the identity in (2.3). □ 

Remark 2.2. It is worthwhile to mention that transforming to the random variable 
I/AtiI-ct) (2.5) is equivalent to the otherwise non-obvious integral identity 

n Ji U/-1) 71 
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This leads to interesting results when the density £, a F x ( x ) has a known form. On the 
other hand, we see that one does not need the explicit density of M a (Fx) to obtain 

the density of M\{Fxy„) = Pa,i~aMa{Fx)- In fact, owing to our goal of yielding simple 
densities for many Dirichlet means from one mean, we see that the effort to calculate the 
density of M\(FxY a ) for each < a < 1 is no more than what is needed to calculate the 
density of Mi (Fx). 

We now see how this translates into the usually difficult problem of explicitly describing 
the density of the finite-dimensional distribution of a subordinator. In the next result, 
we write, for some set C, 

poo 

Ce(C):= l(seC)(e(ds). 
Jo 

Proposition 2.1. Let (&(*);*< 1/6*) denote a GGC(9,F X ) subordinator on [0,1/8]. 
Furthermore, let {C\, . . . ,Ck) denote an arbitrary disjoint partition of the interval [0, 1/9]. 
The finite- dimensional distribution (Ce(Ci), ■ • ■ , Ce(Cfc)) then has a joint density 

k 

Y[9,„F x (x t ), (2.7) 

i=l 

where each o~i = 9\Ci\ > and 2j=i a i = I- The density g ai ,F x * s given by (2.6). That is, 
Ce{Ci) = G\Mi{FxY a . ) and these are independent for i = 1, . . . , k, where M\(FxY a . ) = 
@<ji.i~aiMai(Fx) has density 

-x^- 1 sm(nFxY a (x))c^'' i ' F x^). 

Proof. First, since (Ci, . . . , Ck) partitions the interval [0, 1/9], it follows that their sizes 
satisfy < \C k \ < 1/9 and £ti \C k \ = 1/9. Since Ce is a subordinator, the independence 
of the Ce(Ci) is a consequence of its independent increment property. In fact, these are 
essentially equivalent statements. Hence, we can isolate each Ce(Cj). It follows that for 
each i, the Laplace transform is given by 

E [ e -ACe(Ci)] =e -e\Ci\^ Fx {\) =e -<^i=>(A) ; 

which shows that each Cfl(Ci) is GGC((7i,Fx) for < <Tj < 1. Hence, the result follows 
from Theorem 2.2. □ 

3. Exponential tilting/Esscher transform 

In this section, we describe how the operation of exponential tilting of the density of 
a GGC(#, Fx) random variable leads to a non-trivial relationship between a mean func- 
tional determined by Fx and 9, and an entire family of mean functionals indexed by an 
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arbitrary constant c > 0. Additionally, this will identify a non-obvious relationship be- 
tween two classes of mean functionals. Exponential tilting is merely a convenient phrase 
for the operation of applying an exponential change of measure to a density or more gen- 
eral measure. In mathematical finance and other applications, it is known as an Esscher 
transform and is a key tool for option pricing. We mention that there is much known 
about exponential tilting of infinitely divisible random variables and, in fact, Bondes- 
son [3], Example 3.2.5, explicitly discusses the case of GGC random variables, albeit 
not in the way we shall describe it. In addition, examining the gamma representation 
in (1.1), one can see a relationship to Lo and Wcng [34], Proposition 3.1 (see also Kiichler 
and Sorensen [30] and James [22], Proposition 2.1), for results on exponential tilting of 
Levy processes). However, our focus here is on the properties of related mean functionals, 
which leads to genuinely new insights. 

Before we elaborate on this, we describe generically what we mean by exponential 
tilting. Suppose that T denotes an arbitrary positive random variable with density, say, 
fx- It follows that for each positive c, the random variable cT is well defined and has 
density 

-fr(t/c). 
c 

Exponential tilting refers to the exponential change of measure resulting in a random 
variable, say T c , defined by the density 

_ c-Hl/c)f T (t/c) 
Jf ^>~ E[e" cT ] ■ 

Thus, from the random variable T, one gets a family of random variables (T c ;c > 0). 
Obviously, the density for each T c does not differ much. However, something interesting 
happens when T is a scale mixture of a gamma random variables, that is, T = GgM 
for some random positive random variable M independent of Gg. In that case, one can 
show, see [23], that T c = GgM c , where M c is sufficiently distinct for each value of c. We 
demonstrate this for the case where M — Mg(Fx)- 

First, note that, obviously, cMg(Fx) = Mg(F c x) for each c > 0, which, in itself, is not 
a very interesting transformation. Now, setting Tg = GgMg(Fx) with density denoted 
90.Fx> the corresponding random variable Tg^ c resulting from exponential tilting has 
density 

e- t (l/c) fltf , Fjt (*/c)e°^W (3.1) 



and Laplace transform 

E -S(# x (c(l+A))-^(c)] _ ( 32 ) 



E r e -c(l+A)G e M e (F x ) 



Now, for each c > 0, define the random variable 

d cX 



(cX + 1) 
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That is, the cdf of the random variable A c can be expressed as 



F Ac {y)=Fx 



( 



c(i-y) 



) 



for < y < 1. 



In the next theorem, we will show that M# (fx ) relates to the family of mean functionals 
(Mg(FA c );c > 0) by the tilting operation described above. Moreover, we will describe the 
relationship between their densities. 

Theorem 3.1. Suppose that X has distribution Fx and for each c > 0, A c = cX/ (cX + 1) 
is a random variable with distribution Fa c - For each 6 > 0, let Tg = GgMg(Fx) denote 
a GGC (9, Fx) random variable having density gg,F x - Let Tg^ c denote a random variable 
with density and Laplace transform described by (3.1) and (3.2), respectively. Tg^ c is then 
a GGG{9,Fa c ) random variable and hence representable as GgMg(FA c ). Furthermore, 
the following relationships exist between the densities of the mean functionals Mg(Fx) 
andMg(F Ac ): 

(i) supposing that the density of Mg(Fx), say t;gF X 7 * s known, then the density of 
Mq(Fa c ) is expressible as 



for0<y<l; 

(ii) conversely, if the density of Mg(FA c ), £,9F A * s known, then the density of 
Mg(Fx) is given by 



Proof. We proceed by first examining the Levy exponent [i/jf x ( c (^ + ^)) ~ V'f x (c)] as- 
sociated with Tg. c as described in (3.2). Note that 



and ipF x ( c ) i s °f the same form with A = 0. Hence, isolating the logarithmic terms, we 
can focus on the difference 






log(l + c(l + A).t) - log(l + ex). 



This is equivalent to 
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showing that Tg^ c is GGC(9, Fa c )- This fact can also be deduced from Proposition 3.1 in 
Lo and Weng [34] . The next step is to identify the density of Mg (Fa c ) in terms of the 
density of Mg(Fx)- Using the fact that To = GgMg(Fx), one may write the density of 
Tg in terms of a gamma mixture as 



:(*) 



1 /'OC 



-t/m m -e 



m !;gF x {m)dm. 



m Jo 

Hence, rearranging terms in (3.1), it follows that the density of Tg yC can be written as 



qHf X (c) 



1 poo 



m J 



< cm+ ^^ cm \cm)~ e ^ Fx (m) dm. 



Now, further algebraic manipulation makes this look like a mixture of a gamma(#, 1) 
random variables, as follows, 



t" 



1 poo 



r(0) Jo 



-t(cm+l)/ (cm) 



cm + 1 



(1 + cm) 



dm. 



Hence, it is evident that Mg(FA c ) has the same distribution as a random variable 
cM/{cM + 1), where M has density 



e^' c '(l+H'% fl (m). 
Thus, statements (i) and (ii) follow. 



□ 



3.1. Tilting and beta scaling 

This section describes what happens when one applies the exponential tilting operation 
relative to a mean functional resulting from beta scaling. Recall that the tilting operation 
applied to GgMg(Fx) described in the previous section sets up a relationship between 

Mg(F x ) and Mg(FAc)- Consider the random variable f3g<j,8(i-cr)Mga{Fx) = Mg(FxY„)- 
Then, tilting GgMg{FxY„) as in the previous section leads to the random variable 
GgMg(F c xY a /(cXY„+i)) and hence relates 

Pe*,9(i-*)Me*(F x ) = Mg(F XY J 
to the Dirichlet mean of order 8, 

Mg(F c xYv/(cXY„+l))- 

Now, letting Fa c y„ denote the distribution of A c Y a , one has 
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and hence 



Mg(F cXYJ(cXYa+1) ) 4 Mg{F AcY J = /W(i- CT )Af eCT (F A J. 



(3.3) 



In a way, this shows that the order of beta scaling and tilting can be interchanged. 
We now derive a result for the cases of M\{Fxy c ,) = Pa.i-aM (y {Fx) and Mx{Fa c y^) — 
/3 (T! i_ (J M (T (i 7 ' J 4 c ), related by the tilting operation described above. Combining Theorem 2.2 
with Theorem 3.1 leads to the following result. 

Proposition 3.1. For each < a < 1, the random variables M\(Fxy„) = P<7,i-aM<j{Fx) 
and Mi(FA c Y a ) = /3a,i-a-M rT (FA c ) satisfy the following: 



(i) the density of Mi(Fa c y„) is expressible as 

TCF X Y„ 



e <7l/>F x (c)yCT-l 

(y) = — — — sm 



nc a (l — y) 



<i-y) 



„-<7&F x (v/(c(l-y))) 



forO<y<l; 
(ii) conversely, the density of M\ {F X Y a ) is given by 



£fw (x) 



e -crip Fx (l) x a-l 
7t(l +X) 



l + x 



-CT<S> F (x/(l+x)) 



Proof. For clarity, statement (i) is obtained by first using Theorem 3.1, which gives 

«^(») = ie*'- W (l-ri-W( ;f jL) 

for < y < 1. It then remains to substitute the form of the density (2.5) given in Theo- 
rem 2.2. Statement (ii) proceeds in the same way, using (2.6). □ 



Note that even if one can calculate $ f a for some fixed value of c, it may not be so 
obvious how to calculate it for another value. The previous results allow us to relate their 
calculation to that of $_f x , as described next. 

Proposition 3.2. Set A c = cX/(cX + 1) and define $F Ac (y) =E[log(\y- A c \)I(A c ^y)]. 
Then, for < y < 1, 

®f Ac (y) = $F X ( c ^ V _ y y ) - (c) + log(c(l - y)). 

Proof. The result can be deduced by using Proposition 3.1 in the case a = 1. First, note 
that sin(TCFx( c ri_ y \ )) = shi(nFA c {y)) ■ Now, equating the form of the density of Mi(Fa c ) 
given by (1.5) with the expression given in Proposition 3.1, it follows that 

p -*^ r (y) _ e^*W $ (y/(c(l-y))) 

c(l-y) 
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which yields the result. □ 

Remark 3.1. Wc point out that if G K represents a gamma random variable for n / 6, 
independent of Mg(Fx), then it is not necessarily true that G K M$(Fx) is a GGC random 
variable. For this to be true, Me (Fx) would need to be equivalent in distribution to some 
M k (Fr). In that case, our results above would be applied for a GGC(k,Fr) model. 



4. Examples 

In this section, wc will demonstrate how our results in Sections 2 and 3 can be applied 
to extend results for two random processes recently studied in the literature. The first 
involves a class of GGC subordinators that can be derived from a random mean of a two- 
parameter Poisson-Dirichlct process with a uniform base measure, which was studied 
as a special case in James, Lijoi and Priinster [25]; see Pitman and Yor [37] for more 
details of the two parameter Poisson-Dirichlct distribution. The second involves a class of 
processes recently studied in Bertoin et al. [2]; see also Maejima [36] for some discussion 
of this process. A key component will be the ability to obtain an explicit expression for 
the respective $>f x ■ I n the first example, we do not have much explicit information on 
the relevant density, £sf x ; however, we can rely on a general theorem of James, Lijoi 
and Priinster [25] to obtain . In the second case of the models discussed in Bertoin 
et al. [2], this theorem apparently does not apply. However, we will be able to use an 
explicit form of the density, obtained for a particular value of 9 by Bertoin et al. [2] , to 
obtain $ p x . 

As we shall show, both of these processes arc connected to a random variable Z a , 
whose properties we now describe. For < a < 1, let S a denote a positive a-stable 
random variable specified by its Laplace transform 

E[e" AS <*] =e" AQ . 



In addition, define 

where S' a is independent of S a and has the same distribution. The density of this random 
variable was obtained by Lamperti [31] (see also Chaumont and Yor [4], Exercise 4.2.1) 
and has the remarkably simple form 

sin(7ta) 1 

k (v) = — - — - -T—z 7 — r—7 for V > °- 

Tta y z + 2ycos(7ta) + 1 

Furthermore (see Fujita and Yor [18] and (James [24], Proposition 2.1), it follows that 
the cdf of Z a satisfies, for z > 0, 



376 



L. F. James 



and 



F Za (l/z) = l-±cot~^ COS{na) + 1/z 



not \ sin(7ta) 



1 _■, / cos(7ta) + , 
cot 



Tien \ sin(7to;) 

= l-F Za (z), 

sin (naF Za (z )) = zsm (na 1 - F z „(*))) = — J— ^ (4.1) 

[z 2 + 2z cos(Tta) + l] v " 



sin(27ta) + 2zsin(7ta) 

sin 27ta 1 - F Za (z)]) = — — - 

1 + 2zcos(7Ta) + z l 

2sm(7ta)[cos(7ra) + z] 



1 + 2zcos(7to:) + z 2 

When a = 1/2, 

sm(7t[l-F Zl/2 (z)]): 



(4.2) 



z 2 + l 



4.1. Subordinators derived from an example in James, Lijoi and 
Priinster 



For < a < 1 and 9 > —a, we define the special case of a two-parameter Poisson-Dirichlet 
random probability measures as 

oo k— 1 

^(•)=E^Il( 1 - y ^(-)' 

k=l i=l 

where C/fe are i.i.d. Uniform[0,l] random variables and the are a sequence of inde- 
pendent P a ,B+ka random variables, independent of (Uk)- So, in particular, these random 
variables satisfy E[P Qi g(-)] = Fu(-), where U denotes a Uniform[0, 1] random variable. In 
addition, Po.e is a Dirichlet process. Then, consider the random means given as 

oo fc — 1 „ 1 

M a> e{F u ):=V a> e = Y^U h V h T\{l-V i )= / uP afi {du). 
fe=i i=i - 70 

The \J a ,9 represent a special case of random variables representable as mean functionals 
of the class of two-parameter (a, 0) Poisson-Dirichlet random probability measures - that 
is to say, random variables M a j(Fx), where Fx is a general distribution. An extensive 
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study of this larger class was conducted by James, Lijoi and Priinster [25]. In regards to 
U aj g, they show that U Qj o has an explicit density 

sin(Tta) (a + l)t a (l-t) a 



an [t2a+2 + 2ta+i(l - t) a+1 cos(Tta) + (1 - t) 2a+2 ] ' 
Furthermore, from James, Lijoi and Priinster [25], Theorem 2.1, for 9 > 0, 

U M = M e (F Dat0 ). 

This implies that 

GeV a ^ = GgMe(Fu a0 ) 

are GGC((9, F\j a ). Now, from Vcrshik, Yor and Tslevich [40] (see also James, Lijoi and 
Priinster [25], equation (16)), it follows that 



E[e 



A(a + 1) 



9/a 



(A + l) a+1 -1, 
exp(-0E[log(l + AU a ,o)]), 



where this expression follows from the generalized Sticltjcs transform of order —a of a 
Uniform[0,l] random variable, 

E[(l + At/) Q ]= / (l + \x) a dx^ { ' 



\(a + l) 



A description of the densities of U a ,9 for 6 > —a is available from the results of [25]. 
However, with the exceptions of U Qi i and U ai i_ a , their densities are generally expressed 
in terms of integrals with respect to functions that possibly take on negative values. 
Here, by focusing instead on random variables f3e,i-e^Ja,e for < 8 < 1, we can utilize 
the results in James, Lijoi and Priinster [25] to obtain explicit expressions for their 
densities and the corresponding GGC(6,Fv a ) random variables. 

In particular, we will obtain explicit descriptions for the finite-dimensional distribution 
of a GGC(a, F\j a ), say (Y Q (i),t > 0), subordinator, where T Q (1) = G Q U QjQ and hence 

E[e -AT,(i) ]= *(q + l) , 
L J (A + 1) Q+1 -1 

Although not immediately obvious, one can show that 

z y(oc+i) / / t \ a+l 

U «.o = zl/ a {a+1) - 1 and hence F Va (t) = F Za \^— 

From this, due to the tilting relationship discussed in Section 3, we see that we can also 
obtain results for the GGC(a,F z i/(o+i)) subordinator, say (T* (t),t > 0). To the best of 
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our knowledge, this process and its mean functionals M((F i/(„+i) have not been studied. 
Now, from James, Lijoi and Priinster [25], Theorem 5.2(iii), it follows that 

[ t 2«+2 + 2^+1(1 - t ) a+1 cos(Tta) + (1 - t) 2a+2 f i2a) 

This, combined with our results, leads to an explicit description of the finite- 
dimensional distribution of the relevant subordinators. 

Theorem 4.1. Consider the GGC(a,Fy aa ) subordinator (T a (t),t < l/a) and the 
GGC(a, F z i/( a +i) ) subordinator (T^(i),i < l/a). Let (Ci,...,Cfc) denote an arbitrary 
disjoint partition of the interval (0,1/a] with lengths \Ci\ and set Oi = a\Ci\ for i = 
1, . . . , k. The following results then hold: 

(i) The finite dimensional distribution of (T Q (Ci), . . . , T a (Cfe)) is such that each 
Ya(Ci) is independent and has distribution 

T a (G\) ^GMFy^,), 

where Mi(i 7 y iT .u Q ) = /3 rTu i^ <7i V a . ai . Furthermore, for any fixed < cr < 1, the 
density of Mi(Fy a v a a ) * s given by, for < y < 1, 

(a + l)* / V- 1 BinMl--Pii < „o(v)]) 



[ y 2«+2 + 2^+1(1 - yf +1 COs(TTa) + (1 - y) 2 «+ 2 ] ' 

(ii) The finite- dimensional distribution of (T|,(Ci), . . . , T|,(Cfe)) is suc/i £/iai eac/i 
T|(Ci) is independent and has distribution 

where M\(Fy z i/(o+i)) = f3 ait x- ai M ai (F z i/( Q+ i) ). Furthermore, for any fixed < 
cr < 1, i/ie density of Mi(F y z i/( Q +i)) zs given by, for x > 0, 

a*- 1 (x + i)*(i+°0A»-i s in(7To-[l - F Zq (x a+1 )]) 
[x 2a + 2 + 2x a + 1 cos(Tta) + Ifl^ 

Proof. Statement (i) follows from Theorem 2.2 and Proposition 2.1 in combination 

with (4.3). Noting the relationship between zU^^ and U Qj 0j statement (ii) follows 
from Theorem 3.1(h). □ 



From this, combined with an application of (4.1), we obtain a description for the 
densities of T|(l) and T Q (1). In addition, for a < 1/2, we obtain a description of the 
density of T Q (2) using (4.3). 
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Proposition 4.1. Let T Q (1) and Tjj(l) denote GGG random variables with parameters 
(a,Fjj aQ ) and (a,F z i/( a +i), respectively. Then: 

(i) T Q (1) = GiMi(FY a xs a0 ), where Mi (i*y a u a ) = Pa,i-a^a,a has density, for < 

y<i, 



n [ y 2cc+2 + 2ya +l (1 _ yf + ' cOs(TTa) + (1 - V y a+Z ] 

(ii) T*(l) = G\Mi{F Y z i/(a+i>), w/iere M^i*!^ z i/(o+i)) ftas density 
sin(7Ta) (1 + x) a 



7T [x 2a + 2 + 2x a+1 cos(7ta) + 1] 



/or x > 0. 



(iii) Supposing that a < 1/2, £/ien £/ie GGC(2a, F z \/( a +i) random variable T| a (l) = 
T*(2) is equivalent in distribution to G\Mi(F Y ^ z i/(a+it), where Mi (F y ^ z i/(a+i)) 
has density 

2x 2a - 1 (x + if a+1 sin(7ta)fcos(7to;) + x Q+1 l 

: : 1 il i i i for x > 

[x 2a + 2 + 2x a + l cos(7ta) + l] 2 

4.2. An example connected to Diaconis and Kemperman 

Note that we have the following convergence in distribution results, as a — > : 
M afi {Fu) = M e (F Va ) A Mg(Fu) for > 



and 



u Q , Ai-u. 



Furthermore, setting W = (1 — U)/U = Gi/G[, we have 

M e (F z ./(«+i)) AM e (F w ) and -4 W, 

where the last statement can be read from Chaumont and Yor [4], page 155 and page 
169. It is then natural to investigate the laws of the random processes connected with 
the GGC(#, Fjj) and GGC (6*, Fw) random variables. It is known from Diaconis and 
Kemperman [9] that the density of M\(Fu) is 

-sm(TCy)y~ v {l -y)~ (1 ~ v) for0<y<l. (4.4) 

7T 

Note, furthermore, that T\ = GiM\{Fjj) is GGC(1, Fjj) and has Laplace transform 



E [ e -AG 1 M 1 (F [ ;)] =e -^(A) =e (l + A) 



-(0+i)A) 
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Now, GiMi(F w ) is a GGC(l,F w ) with ipF w (X) = ^-log(A). Theorem 3.1 shows that 
Mi(Fjj) arises from tilting the density of GiMi(F\y)- The density of Mi(Fw) is obtained 
by applying statement (ii) of Theorem 3.1 to (4.4), or by statement (ii) of Proposition 
3.1, and is given by 



€f w (x) = - sin 

7T 1 



TlX 



-x/(l+x) 



for x > 0. 



We now apply Theorem 2.2 and Proposition 2.1 to give a description of the finite- 
dimensional distribution of the subordinators associated with the two random variables 
above. 



Proposition 4.2. Let U denote a uniform [0,1] random variable and let W = Gi/G[ 
denote a ratio of independent exponential (1) random variables. 

(i) Suppose that (d(t);0 < t < 1) is a GGC(1, Fjj) subordinator. Then, for (Ci, ... ,Ck), 
a disjoint partition of (0, 1), the finite- dimensional distribution has joint density 
as in (2.7), with 

Q.uFu fo) = / e-**/"— sin(7ra,(l - yM^-^l - y)-^ 1 ^ dy 
Jo n 

for i=l,...,k. 

(ii) That is, Ci (Ci) = GiMi^Fuy^ .) and they are independent for i = 1, ...,k. Fur- 
thermore, the density of M\(FuY a .) = P, 7i _i- ai M c , i {Fii) is 

— sin(7t^(l - yW^-^il - j,)"^ 1 -") 

7T 

/or < y < 1. 

(iii) // (Ci(t); < t < 1) is a GGC(l,F\y) subordinator, then the finite- dimensional 
distribution (£i(Cx), . . . , Ci(Cfc)) is noio described, with 

g ai iFw (*,) = r e-*«/™ - Bin f W^)" 2 d™. 

(iv) TTiai is, Ci(Ci) = GiAfi(i i V"K (r . ) an-rf i/ie?/ are independent for i = 1, . . . ,k. Fur- 
thermore, the density of Mi(F\yY„.) = Pa-i ,1-0-4 M^iFy/) is 

7T \1 +Xj 

for x > 0. 
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Proof. This now follows from Theorem 2.2, Proposition 2.1 and (4.4). Specifically, note 
that C(Fij) = (0, oo), so for any < a < 1, 

sin(7ti 7 i7Y" ii («)) = sin(7tcr(l — u)) 

for < u < 1 and otherwise. Furthermore, from (4.4), or by direct argument, it is easy 
to see that 

<l> Fu (y) = -log(y-y(l-y)- {1 - y) )-l. 

This fact is also evident from Diaconis and Kemperman [9]. It follows that Mi(Fjjy„) 
has density 

— sin(7ta(l - y))y CT ( 1 ^- 1 (l - yy a(1 ' y) for < y < 1. 
The density for Mi(Fwy„) is obtained in a similar fashion by Proposition 3.1. □ 
Remark 4-1- Setting 

. d. cGi 



cG x + G[ 



one can easily obtain the density of the random variable Mi (Fa c ) for each c > by us- 
ing statement (ii) of Theorem 3.1. Note, also, that one can deduce from the density of 
Mi(F\y) that $f w (x) = [x/(l + x)] log(a;). Hence, in this case, an application of Propo- 
sition 3.2 shows that 

c(l-y)+y \c(l-y)J c-1 
We note that, otherwise, it is not easy to calculate $^ c , in this case, by direct arguments. 



4.3. The finite-dimensional distribution of subordinators of 
Bertoin et al. 

Our final example shows how one can apply the results in Sections 2 and 3 to obtain 
new results for subordinators recently studied by Bertoin et al. [2]. In particular, they 
investigate properties of the random variables corresponding to the lengths of excursions 
of Bessel processes straddling an independent exponential time, which can be expressed 
as 

w(«) _ ( Q ) 

where, for any t > 0, 

g { t a) :=s\xp{s<t,R s = Q}, d t (a) := inf{s > t,R s = 0} (4.5) 
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for (Rt,t > 0) a Bessel process starting from with dimension d = 2(1 — a), with < d < 2 

or, equivalently, < a < 1. Additionally, e = G\, an exponentially distributed random 
variable with mean 1. See also Fujita and Yor [19] for closely related work. 

In order to avoid confusion, we will now denote relevant random variables appearing 
originally as A a and G a in Bcrtoin et al. [2] as E Q and G Q , respectively. From Bertoin 
et al. [2], let (£ Q (i);t > 0) denote a subordinator such that 

E [ e -AE„(t)] = ((A + 1) Q - A a )* 

= exp(-t(l - a)Epog(l + A/G a )]), 

where, from Bertoin et al. [2], Theorems 1.1 and 1.3, G a denotes a random variable such 
that 



1 + zl'_ a c 



and has density on (0,1) given by 

asin(7ta) u a_1 (l — u) c 



(1 - a)7T u 2a _ 2(1 _ COs(7ta) + (1 - u) 



2a ■ 



Hence, it follows that the random variable 1/G Q takes its values on (l,oo) with proba- 
bility one and has cdf satisfying 

1 - F 1/Ga (x) = F Ga (l/x) = F Zl _ a ((x - iy a ). 

As noted by Bertoin et al. [2], (E a (t);t > 0) is a GGC(1 — a,-Fi/G a ) subordinator, where 
the GGC(1 — a,Fi/ Ga ) random variable E Q = E Q (1) satisfies 



a ) A Gi-q, Gi_ c 



where t/ denotes a uniform[0, 1] random variable and, for clarity, G\~ a is a gamma(l 
a, 1) random variable. This means that the density of E Q is 

-ar a-1 (l-e- x ) forrOO. 



r(i-a) 

It is evident, as investigated in Fujita and Yor [18], that 

Ml _ a (F 1/G jiJ-^C7-V« 

Pa.l 
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Remark 4-2. Note that when a = 1/2, G1/2 = /?i/2,i/2- It is known that for each fixed 
t, 



Si/a(t) 



Gt/2 



/3l/2,(l+t)/2 ' 

which implies that 

M t/2 (F 1/Gl/2 )=M t/2 (F 1/Pl/21/2 ) 4 1 . (4.6) 

Pl/2,(l+f)/2 

This result may be found in James and Yor [27]. Related to this fact, Cifarclli and 
Melilli [7] have shown that M t / 2 (F^ 1/21/2 ) = /3(t+i)/ 2 ,( t +i)/2 for i > 0. 

In regards to exponentially tilting GGC(1 — a,Fi/ Ga ), note that for c > 0, 

c/G„ c 



Thus, a GGC(1 — a, F c /(<s a + c )) subordinator, say (E* c (i),t < 1/(1 — a)), arises from 
exponential tilting. Naturally, the density of E^ c (l)/c is given by 



[(c+l) Q -c Q ]r(l-a) 



Equivalently, E|, jC (l) = Gi_ Q Afi_ a (F c/(Ga+c )), where Mi_ a (F c /( Gct+c )) has density 



-u for < u < 1. 



(c+l) Q -c a c+1 

Now, using the facts discussed above, we will show how to use the results in Section 2 to 
explicitly describe the finite-dimensional distribution of the subordinators (E Q (f),t > 0) 
and (E* c (t),t > 0) over the range <t < 1/(1 — a). Additionally, the analysis will also 
yield expressions for mean functionals based on F\/ Ga . First, note that, using (2.4), one 
has 

. , , Jsin(7t(l -a)F Ga (l/x)), ifx>l, 

sm(TtF Yl _ a/Gct (x )) = i . (4.7) 
L sin(7t(l — a)), if < x < 1, 

where, again using the properties of Fz 1 _ a , as deduced from James [24], Proposition 
2.1(iii), 

, , , , ,, sin (nil — a)) 
sm(7i(l-a)F Ga (l/x)) = U U (4.8) 

[(x - 1) 2Q - 2(.x - 1) Q cos(7ta) + 1] 

We now use this to calculate 

'!>/ . (a?) = E[log(|x - l/G a \)l{x + 1/G a )]. (4.9) 
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Proposition 4.3. For < a < 1, consider $>F 1/Ga (x) as defined in (4-9). Then, 

/- 1 r 

log 



[(a; - 1) 2q - 2(x - 1) Q cos(Tta) + 1] 



(4.10) 



ifx>l, 
if0<x< 1. 



= J 2(1 -a) 

J-iog^/ii-a-^n), 

1 — a 

Proof. Using simple beta-gamma algebra, we have 

Hence, applying Theorem 2.1, with = 1 and a = 1 — a, it follows that S Q is also 
GGC(l,F Yl _ / G J and 



By standard calculations, the density of £> Q = Pi- a ,a/ Pa,i is given by 

sin(7t(l-a)) _ x 



(4.11) 



71 



However, we see from (4.11) that B a = Mi(F yi _ Q / Ga 
the density of B a can be written as 

x~ a 

fB a (x) = sm(nF Yl _ a/Ga (x))e 



x-v-^l-il-xFlixKl)]. 

). Hence, Theorem 2.2 applies and 



-(l-a)*p 1/G (a) 



Now, equating the two forms of the density of B a and using (4.7) and (4.8), one then 



obtains the expression for <&p 1 , G . 
Now, for z > 0, define the function 



□ 



S a . a {z) = sin(7tcr J F Zl _ Q {z~ a ))[z Za - 2z a cos(Tta) + 1] 



r/(2(l-c*)) 



and define, 



I l-V.cr^ — 1): 



u-i<7 / (1 — a) 



\ix<l, 
if x > 1. 



Hereafter, (C±, . . . ,Ck) will denote an arbitrary disjoint partition of the interval 
(0, 1/(1 — a)] with lengths |C,|, and u, = (1 — a)|Ci| for i = 1, . . . , fc. 
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Theorem 4.2. Consider the GGC{\ — a, -F\/g q ) subordinator (T, a (t),t < 1/(1 — a)) and, 
for each fixed c> 0, the GGC(l — a,F c /(c a+c ) subordinator (Ej, c (t),t < 1/(1 — a)). The 
following results then hold: 

(i) The finite- dimensional distribution of (£ Q (Ci), . . . , £ a (Cfc)) is such that each 
S Q (Ci) is independent and has distribution 

E a (Ci)^(3iMi(F ni/G J, 

where Mx(-Fy„./G a ) = /^.l-o-j-Ma-^-Fi/Ga)- Furthermore, for any fixed < o" < 1, 
i/ie density of Mi(Fy a /& a ) is given by 



l x -(<ra/(l- a )+l) Va M forx>Q. 
71 

iii) For the GGC(1 — a,F c /rig a+c )) process, SJ, c; if follows that each 

4, c (a) = GWi(^ iC/(Ga+c) ), 
where for each < <r < 1, Mi(FY aC /(G a+c )) has density 
[c{l-y)r' {1 ~^V a (y/{c{l-y))) 



7t[(c + 1) Q - c Q ]' T y CTQ /(i-«)+i 



/or < y < 1. 



Proof. From Theorem 2.2, we have that the general form of the density of Mi(F Yt7 /G a ) 
is given by 

— sin(7tF^ /G >))e-^/ Gc >). 
The proof is completed by applying Proposition 4.3 and (4.7) and (4.8). □ 

Remark 4-3. The process E a c (t)/c is well defined for c > and < a < 1, and presents 
itself as an interesting class worthy of further investigation. Letting c — > 0, it is evident 
that ££, c (l)/c converges to S Q (1). As shown by Bcrtoin et al. [2], Section 3.6.3, £o,c(l)/c, 
for c > 0, has a similar interpretation as S a (l), but where the Bessel process (Rt,t > 0) 
is now replaced by a diffusion process whose inverse local time at is distributed as a 
gamma subordinator (ji/c; I > 0). Furthermore, albeit not explicitly addressed in Bcrtoin 

et al. [2], the random variable £ Q . c (l)/c = di Q ' c ' — ge has a similar interpretation 
where (Rt , t > 0) is now replaced by a process (R { t a ' c \t> 0) whose inverse local time 
is distributed as a generalized gamma subordinator, that is, a subordinator whose Levy 
density is specified by Cy~ a ~ 1 e~ cy for y > 0. This interpretation may be deduced from 
Donati-Martin and Yor ([10], see page 880 (l.c)), where i?( Q ' c ) equates with a downwards 
Bessel process with drift c. 
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Bertoin et al. [2] also show that a GGC (1 — a,Fc a ) random variable satisfies 

Gi- a Mi- a {F<Q a ) = G\- a U. 

Hence, the Laplace transform of the GGC (1 — a, F& a ) subordinator, say (Z^ 1 {t),t>0), 
is given by 



Additionally, using the fact that 



J-4 * +i^U 1/ "+i 



leads to 



Z l-a 



M r - a {F z y a ) = M 1 _ a {F yGa ) - 1 = ^-, 



which leads to a description of a GGC(1 — a, F 7 i/ a ) subordinator. The above points may 
also be found in the survey paper of James, Roynette and Yor [26]. 

Theorem 4.3. Consider the GGC(l — a, F i/ a ) subordinator (Z a (t),t < 1/ (1 — a)) and 

the GGC(l — a,F(; a ) subordinator (Z^ 1 (t),t < 1/(1 — a)). The following results then 
hold: 

(i) The finite- dimensional distribution of {Z a {C\) 1 . . . , Z a (Ck)) is such that each 
Z a (Ci) is independent and is equivalent in distribution to 

Z a (C i ) = G 1 M 1 (F YirtZ i /J .J. 
Furthermore, for any fixedO < a < 1, the density of M\{F _i/ a ) = j3 a _i-aM a {F _i/ a ) 

17 1 — a. X — ot 

is given by, for z > 0, 



n{l + zY ni - a) 

achZ\ L (C i )=G 1 

has density 



(ii) Similarly, each (C;) = GiAfi(i 7 y <J G a ) an d, for each fixed < a < 1, Mi(i 7 V (T G Q ) 



2 ir- 1 (i-y) ra/<1 - a) 5 Q , <r ( T 



.'7 



Proof. Apply Theorem 2.2 and Theorem 3.1, where, from (4.12), 
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Remark 4-4- Note that as a — > 1, 

M 9 (F c jAM4f[/) and Mj(^/. ) 4 M fl (.FW). 

Hence, they have the same limiting behavior, described in Section 4.2, as the random 
variables in Section 4.1. 
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